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■(f^ \ Abstract. We have developed a physical model for the non-perturbative x-shape of parton density 

f*^) , functions in the proton, based on Gaussian fluctuations in momenta, and quantum fluctuations of 

the proton into meson-baryon pairs. The model describes the proton structure function and gives 
Cn| \ a natural explanation of observed quark asymmetries, such as the difference between the anti-up 

and anti-down sea quark distributions and between the up and down valence distributions. We also 
find an asymmetry in the momentum distribution of strange and anti-strange quarks in the nucleon, 
large enough to reduce the NuTeV anomaly to a level which does not give a significant indication 
of physics beyond the standard model. 
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The low-scale parton density functions give a description of the hadron at a non- 
perturbative level. The conventional approach to these functions is to make parameter- 
izations using some more or less arbitrary functional forms, based on data from deep 
inelastic scattering and hadron collision experiments. Another approach, however, is to 
start from some ideas of the behavior of partons in the non-perturbative hadron, and 
build a model based on that behavior. The advantage with this approach is that the suc- 
cesses and failures of such a model allows us to get insight into the non-perturbative 
QCD dynamics. The model presented here, and described in detail in flyD, describes 
the F2 structure function of the proton, as well as sea quark asymmetries of the nucleon. 
Most noteworthy, our model predicts an asymmetry between the momentum distribu- 
tions of strange and anti- strange quarks in the nucleon of the same order as the newly 
reported results from NuTeV Jj]. 

This work extends the model previously presented in The model gives the four- 
momentum k of a single probed valence parton (see Fig. [T^ for definitions of mo- 
menta) by assuming that, in the nucleon rest frame, the shape of the momentum dis- 
tribution for a parton of type i and mass m,- can be taken as a Gaussian fj(k) = 
N(Oi, mi) exp { — [(ko — mi) 2 + k 2 + k 2 + k 2 ] /2a 2 } which may be motivated as a result 
of the many interactions binding the parton in the nucleon. The width of the distribution 
should be of order hundred MeV from the Heisenberg uncertainty relation applied to the 
nucleon size, i.e. <7, = l/d^. The momentum fraction x of the parton is then defined as 
the light-cone fraction x = k + / p + . We impose constraints on the final-state momenta in 
order to obtain a kinematically allowed final state, which also ensures that < x < 1. Us- 
ing a Monte Carlo method these parton distributions are integrated numerically without 
approximations. 
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FIGURE 1. Illustration of the processes (a) probing a valence parton in the proton and (b) a sea parton 
in a hadronic fluctuation (letters are four-momenta), (c) shows the resulting parton distributions at the 
starting scale Qq. 

To describe the dynamics of the sea partons, we note that the appropriate basis for the 
non-perturbative dynamics of the bound state nucleon should be hadronic. Therefore we 
consider hadronic fluctuations, for the proton 

\p) = a \po) + a p7t o\pn ) + a n7t +\nn + } + ... + a AK \AK + } + ... (1) 

Probing a parton i in a hadron H of a baryon-meson fluctuation \BM) (see Fig. Eh) gives 
a sea parton with light-cone fraction x = xnxi of the target proton. The momentum of 
the probed hadron is given by a similar Gaussian, but with a separate width parameter 
Oh- Also here, kinematic constraints ensure that we get a physically allowed final state. 
The procedure gives xh ~ Mr /(Mb + Mm), i-e. the heavier baryon gets a harder spec- 
trum than the lighter meson. The normalization of the sea distributions is given by the 
amplitude coefficients a\ M of Eq. These cannot be calculated from first principles 
in QCD and are therefore taken as free parameters to be fitted using experimental data. 

The resulting valence and sea parton ^-distributions apply at a low scale Qq, and the 
distributions at higher Q 2 are obtained using perturbative QCD evolution at next-to- 
leading order. 

The model has in total four shape parameters and three normalization parameters, 
plus the starting scale, to determine the parton densities u, d, g, u, d, s, s. These are 
(with values resulting from fits to experimental data as described below): 

a u = 230 MeV a d = 170 MeV a g = 77 MeV o H = 100 MeV 

a 2 o = 0.45 a 2 . = 0. 14 a 2 = 0.05 Q = 0.75 GeV (2) 

The resulting parton densities are shown in Fig.dc). 

In order to fix the values of the model parameters, we make a global fit using sev- 
eral experimental data sets: Fixed-target F2 data to fix large-* (valence) distributions 
(Fig.EK); HERA F2 data for the gluon distribution width and the starting scale Qq; d/u- 
asymmetry data for the normalizations of the \pit ) and \mt + ) fluctuations (Fig. 3); 
and strange sea data to fix the normalization of fluctuations including strange quarks 
(Fig. 4a). We have also compared with W ± charge asymmetry data as a cross-check on 
the ratio of Gaussian widths for the u and d valence quark distributions (Fig. Eb)- It is 
interesting to note that this simple model can describe such a wealth of different data 
with just one or two parameters per data set. 
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FIGURE 2. Left: The proton structure function F 2 {x, Q 2 ) for large x values; NMC and BCDMS data 
01 01 compared to our model, also showing the results of ±20% variations of the width parameters a u 
and Gd for the u and d valence distributions. Right: The charge asymmetry for leptons from W ± -decays 
in pp collisions at the Tevatron |7] compared to our model, with best-fit parameters and a 20% reduced 
width of the valence d quark distribution. 

In our model, the shape difference between the valence u and d distributions in the 
proton, apparent from the W ± charge asymmetry data, is described as different Gaussian 
widths. This would correspond to a larger effective volume in the proton for d quarks 
than for u quarks, an effect which could conceivably be explained by Pauli blocking of 
the u quarks. 

Since the proton can fluctuate to n° and n + by \pn°) and \nn + ), but to k~ only by 
the heavier |A ++ 7T~), we get an excess of J over u in the proton sea. Interestingly, the 
fit to data improves when we use a larger effective pion mass of 400 MeV (see Fig. 3). 
This might indicate that we have a surprisingly large coupling to heavier p mesons, or 
that one should use a more generic meson mass rather than the very light pion. 

The lightest strange fluctuation is \AK + ). If we let this implicitly include also heavier 
strange meson-baryon fluctuations, we can fit the normalization ajj, to strange sea 



data (see Fig. 4a). The result corresponds to Jq(xs + xs)dx/ f (xu + xd)dx ~ 0.5, in 



agreement with standard parton density parameterizations. We note that this indicates a 

,2 
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normalization oc 1 /AM B m = 1 / {M B +M M - M p ) rather than the expected oc 1 j hM l B 



Since the s quark is in the heavier baryon A and the s quark is in the lighter meson K + , 
we get a non-zero asymmetry S~ = J dx[xs(x) —xs(x)] in the momentum distribution 
of the strange sea, as seen in Fig. 4b and 5. Depending on details of the model, we get 
the range 0.0010 < 5 < 0.0023 for this asymmetry. 

This is especially interesting in connection to the NuTeV anomaly NuTeV found, 



based on the observable R 
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FIGURE 3. Comparison between our 
model and data from the E866/NuSea col- 
laboration 1 8]: (a) ii(x)/d(x) (b) xd(x) — 
xu(x). The full line uses the physical pion 
mass, while the dashed line uses an effec- 
tive pions mass m eS = 400 MeV as dis- 
cussed in the text. 



X 

+ 

x 



0.4 r 



F CCFR 
L — s/u=1.0 

: — S/U=0.5 
z - ~ s/u=0.0 

| Q 2 =4 GeV 2 


^Q^OGe^ 


[q 2 =100 GeV 2 


-| i t i 1 ? i "TrT^Hrr-fr. 







[T0.015E 




0.1 0.2 0.1 0.2 0.1 0.2 

xxx 



-0.01 5 E 



0.1 0.2 0.3 0.4 0.5 



FIGURE 4. Left: CCFR data 1 10] on the strange sea distribution (xs(x) +xs(x))/2 compared to our 
model based on \AK) fluctuations with different normalizations. Right: The strange sea asymmetry 
s~ (x) = xs{x) — xs(x) (at Q 2 = 20 GeV 2 ) from the model and combined with the function F(x) accounting 
for NuTeV's analysis giving Asin 2 9 W = J dxs~ (x)F(x) = —0.0017. The uncertainty bands correspond 



to the uncertainties for S 
0.05 r— 



and Asin 9w quoted in the text. 




FIGURE 5. Comparison between the strange and charm 
sea obtained from our model with the inclusion of the 
AcD fluctuation. The normalization is here taken to be 
oc \j (M B +Mm —M p ) as suggested by strange sea data. 

deviation of sin 2 9 W compared to the Standard Model fit: sin 2 0$ uTeV = 0.2277 ±0.0016 
compared to sin 2 0^, M = 0.2227 ±0.0004. However, an asymmetric strange sea would 
change their result, since v only have charged current interactions with s and v with 
s. Using the folding function provided by NuTeV to account for their analysis, the s- 
s asymmetry from our model gives a shift —0.0024 < Asm 2 6w = Jq dxs~(x)F(x) < 
—0.00097, i.e. the discrepancy with the Standard Model result is reduced to between 
1.6a and 2.4a, leaving no strong hint of physics beyond the Standard Model. 

We have also considered charmed fluctuations. The lightest charmed baryon-meson 
fluctuation \AcD) gives c and c distributions as in Fig. 5, where the normalization ce 2 ^ 
is taken to be \/AMbm, as suggested by the strange sea normalization. However, in 
order to conform to the EMC F 2 C data at large x, a normalization close to 1 / AM 2 — seems 

to be enough 111 ill . 
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